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A full strength Coulomb interaction between trapped electrons can be felt only in absence of a 
neutralizing background. In order to study quantum degenerate electrons without such a back¬ 
ground, an external trap is needed to compensate for the strong electronic repulsion. As a basic 
model for such a system, we study a trapped electron pair in a harmonic trap with an explicit 
inclusion of its Coulomb interaction. We find the eigenenergy of the ground state, confirming earlier 
work in the context of harmonium. We extend this to a complete set of properly scaled energies 
for any value of the trapping strength, including the excited states. The problem is solved either 
numerically or by making harmonic approximations to the potential. As function of the trapping 
strength a crossover can be made from the strongly to the weakly-coupled regime, and we show that 
in both regimes perturbative methods based on a pair-wise electron description would be effective 
for a many-particfe trapped electron system, which resembles a Wigner crystal in the ground state 
of the strongly coupled limit. 

PACS numbers: 02.50.Ga, 32.80.Rm 


INTRODUCTION 

The strong Coulomb interaction between electrons is 
rarely felt at full strength, as their repulsion is commonly 
screened by the presence of positive charges. This is usual 
fashion in condensed matter physics, where the valence 
electrons can be treated as a weakly-interacting electron 
gas. Such systems are said to be in the weakly to in¬ 
termediate coupling regime, which has been successfully 
described by Landau Fermi Liquid Theory. This weakly 
interacting electron gas is also quantum degenerate in 
most metals at room temperature. 

The coupling parameter of these many-particle systems 
is generally considered as the ratio of potential energy 
over kinetic energy. Also plasma physics systems are gen¬ 
erally weakly coupled: here screening is governed by the 
Debye length, which make the effective electron inter¬ 
actions of short range nature. Only in charged particle 
accelerators, the pure Coulomb interaction can be appre¬ 
ciated at full strength. However at the same time, this 
strong repulsion makes it very difficult for an electron gas 
to become strongly-coupled as a result of heating, which 
effectively prevents quantum degeneracy as well. 

Nevertheless, there is an ever higher demand for bright¬ 
ness in the held of charged particle accelerators. Here 
the holy grail would be to achieve ultimate conditions 
of brightness, by creating quantum degenerate electron 
or ion beams. Particle accelerators are heavily used in 
medical and industrial applications and for driving high¬ 
brightness X-ray light sources [1]. Ultimately this would 
require electron beam phase space densities close to de¬ 
generacy, three to four orders of magnitude higher than 
in present state-of-the-art electron microscopy beams [2] . 
In this paper, we investigate some basic properties of 
trapped electrons, which are well-behaved systems where 
heating effects can be controlled, and we investigate dif¬ 


ferent regimes of coupling and degeneracy. The trapping 
frequency uj acts as a control parameter, which allows us 
to study a crossover between different limiting coupling 
regimes. One limit is a strongly-coupled freely propa¬ 
gating electron gas, which could be an electron beam for 
instance. In this regime a Wigner crystal is known to 
exist as the ground state [3]. The other limit, for large 
trapping frequencies, is then a weakly-coupled trapped 
electron gas. 

Our approach allows for new opportunities in pertur- 
batively describing few- and many-particle Coulomb sys¬ 
tems. We offer a simple description of the crossover in 
terms of a reference frequency, that depends on the elec¬ 
tron mass. We extend to other types of charged particles 
including different masses, and discuss the feasibility to 
investigate the crossover in these systems, based on mass 
and trapping possibilities. 

Cold atomic gases offer a promising route towards ulti¬ 
mate brightness conditions. In particular Rydberg crys¬ 
tals are very close to the quantum plasma regime, and 
extraction of the electrons could be done with a mini¬ 
mum amount of heating, such that the electrons could 
be trapped in Paul traps [4], Penning traps [5, 6], and 
Ponderomotive traps [7]. Also other ways towards de¬ 
generacy are currently being explored, such as hybrid 
field electron extraction using femtosecond laser pulses 
[8], direct cooling in ion traps, and sympathetic cooling 
with atoms [9, 10]. 

The quantum-mechanical problem of particles with 
short-range interactions in a harmonic oscillator has been 
well studied [11, 12]. Moreover, appropriate studies 
where the interaction is of the Coulomb type have been 
done as well [13-15]. Identical particles in such a trap 
are referred to as harmonium [16-20], which has been 
studied as an exactly solvable model for the artihcial he¬ 
lium atom, where the interaction of the electrons with 
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the nucleus has been replaced by a harmonic confine¬ 
ment. Analytical solutions have been found for a certain 
set of oscillator frequencies [14] , numerical solutions that 
interpolate between them [17], as well as an analytical in¬ 
terpolation model [18]. Here we demonstrate another set 
of analytical interpolating solutions, and we show that 
the effective Coulomb interaction between the particles 
can in certain cases be considered as a short-range inter¬ 
action. 

The problem is approached by writing the wave func¬ 
tion in the angular momentum basis. In the problem 
of Coulomb scattering this basis is not convenient since 
the phase shift does not converge, in other words, the 
Coulomb scattering problem is never asymptotically free 
[21]. In our problem, however, the harmonic trap cures 
this problem; the eigenstates of the harmonic oscillator 
take over the role of the asymptotic free states. There¬ 
fore, the angular momentum basis is indeed convenient. 
As a result, the eigenstates of the harmonic oscillator 
define our reference energy which will be taken as tiuj. 

PROBLEM DESCRIPTION 


by Eq. (3). As mentioned before, the angular momen¬ 
tum basis is a good basis since the harmonic oscillator 
states serve as the asymptotic free states, which do not 
exist for a simple Coulomb scattering problem. 


EXACT SOLUTIONS 

In general, the Schrodinger equation Eq. (3) is not ana¬ 
lytically solvable, but rather the problem is quasi-exactly 
solvable [23, 24]. This problem has been studied by Taut 
[14], amongst others. Taut has shown that exact solu¬ 
tions can be found by assuming that the wave functions 
can be written as 

u{r) = P{p) exp > (5) 

where p = and 

P{p) = (6) 

i=0 


When considering two trapped electrons, it is conve¬ 
nient to transform to the center of mass and relative co¬ 
ordinates. With this transformation, the problem in the 
center of mass reduces to a three dimensional harmonic 
oscillator, which is easily solved [22], whereas the mo¬ 
tion in the relative coordinate r is governed by a three 
dimensional Schrodinger equation 

H\^)=(^^^+V{r)^\^)=E\^) ( 1 ) 

where p, = ^nie is the effective mass and the potential 


V{r) 


^ 2 2 
-puj r 


Aneor 


( 2 ) 


only depends on the length of the relative coordinate r = 
|r|. ct) is the frequency of the harmonic oscillator and the 
second term is the Coulomb interaction for two electrons. 
Since Eq. (2) is a central potential we can reduce this 
further to a one-dimensional Schrodinger equation which 
reads 


(Pu 
2p dP 


^ 2 2 
-pu! r 


Aneor 


2pP ) 


u = Eu, 


( 3 ) 

by writing the wave function in the angular momentum 
basis 


(r|») = 

r 


( 4 ) 


with quantum numbers I and m as usual. Thus the prob¬ 
lem is reduced from solving a six-dimensional eigenprob- 
lem to solving the one-dimensional eigenproblem given 


is a polynomial of finite degree. With this Ansatz, Eq. (3) 
yields a recurrence relation for the coefficients a^, which 
can be solved when j is given. The corresponding solu¬ 
tions are shown in Fig. 1, where we show the eigenen- 
ergies normalized with respect to the oscillator energy 
hu) versus the oscillator frequency oj. Discussion of the 
eigenenergies is delayed to the next section where full so¬ 
lutions can be found. Taut’s solutions correspond only 
to a countable number of solutions (namely integer j) 
where no solution can be found above a frequency of 
32 ™ 2 °g 2^3 = which is, as we will see below, half the 
Hartree frequency. As we seek solutions for any set of 
parameters we thus need other methods. 


FULL SOLUTIONS 

A complete solution of the problem, i.e. for every set 
of parameters, can be found by solving Eq. (3) numeri¬ 
cally. The eigenenergies are shown as solid lines in Fig. 1 
and match Taut’s solutions as expected. We observe two 
distinct regimes. First, for large frequencies we see that 
the energies are approximately that of a three dimen¬ 
sional harmonic oscillator, or in the limit of w —)■ oo these 
seem to match exactly. We thus can conclude that in 
the regime of large frequencies the Coulomb repulsion is 
ineffective and we call this the free harmonic limit. Ac¬ 
cordingly, for large frequencies the Coulomb interaction 
can be considered for this system as a short range in¬ 
teraction. This is a counter-intuitive result, as a large 
frequency corresponds to a small separation distance r, 
and therefore to a large Coulomb interaction energy. 
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Figure 1: Exact eigenenergies as found by Taut’s method 
(dots) as well as the numerical solutions (lines) of the dif¬ 
ferential equation Eq. (3) as function of the frequency uj. ujq 
is a typical frequency of the problem and is exactly two times 
the highest frequency for which Taut’s method can find an 
exact solution. 

On the other hand, when the frequency is small we see 
that the energy diverges from that of a harmonic oscilla¬ 
tor. Using the definition of the weakly and strongly cou¬ 
pled regimes we can identify this regime as the strongly 
coupled / strongly interacting regime. In this regime the 
Coulomb repulsion is relatively strong. We can relate 
both frequency regimes to densities since a smaller fre¬ 
quency corresponds to a smaller separation distance r 
and thus a lower density. So concluding, for a low den¬ 
sity the interaction is strong and for a high density the 
interaction is weak. This is exactly opposite to what one 
would expect for a short-range interaction force in an 
external field. 


HARMONIC APPROXIMATIONS 


In this section, we derive analytical solutions for this 
problem by using an approximate expression for the ef¬ 
fective radial potential 


VMr) = 


hH{l + 1 ) 


(7) 


This was first done by Taut [14] and in the context of 
harmonium [16, 17] for a partial list of energies, however, 
here we extend this to a complete set of properly scaled 
energies, including the excited states. 

In the free harmonic limit we use 


to find 


U(r) = 

(8) 

E , 3 

(9) 

— =2n+ +-, 


which indeed corresponds to the 3D harmonic oscillator 
energies found for large uj as mentioned in the previous 
section. For small frequencies, thus in the strongly in¬ 
teracting regime, we make a Taylor approximation of the 
potential around its minimum: 


V(r) = Vo + (r - rof . (10) 


For Z = 0 we find the eigenenergies 


hio 2'^/^ Wo/ 




( 11 ) 
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where wq = = Eh/^ is the Hartree frequency, 

which is the typical frequency for this problem. It corre¬ 
sponds to the situation where the minimum of Eq. (7) 
is exactly \/2 times the Bohr radius, ro = v72oo = 
^ • "l'^® term corresponds to the divergence 

of ^ for small frequencies given by a typical exponent - 
1/3 and a proportionality constant directly related to ujq. 
For non-zero I using the same approximation we recover 
the same asymptotic behavior for small frequencies. 

The second term in Eq. (11) corresponds to the spacing 
between different levels. The approximate eigenenergies 
given by Eq. (9) and Eq. (11) can be compared to the 
numerical eigenenergies, see Fig. 2. From Fig. 2 we see 
that the harmonic approximations work well in both the 
strongly and weakly interacting regime. 



Figure 2: Comparison of numerical eigenenergies (solid) with 
the eigenenergies Eq. (9) from the free harmonic approxima¬ 
tions (dash-dot) and the eigenenergies Eq. (11) from the har¬ 
monic approximation around the minimum (dashed). 
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HELLMANN-FEYNMAN 

Using the Hellmann-Feynman theorem [25, 26], given 
here as 


function of the frequency for general quantum numbers n 
and 1. The found energies are presented in Fig. 3, where 
the asymptotic behavior for small frequencies given by 
the first term of Eq. (11) is subtracted for clarity. 


^ f \u{r)f dr = ( 12 ) 

du! ,1 


we are able to calculate when we know the solution 
M(r), which gives us additional information about the 
eigenenergies we are looking for. We assume u{r) is a flat 
wave function on the classically allowed domain of the 
harmonic oscillator potential. This is justified at large 
relative energies and a sufficiently low frequency, which 
corresponds to the strongly interacting regime. With this 
approximation we can determine (r^) for n = I = 0 and 
find 


(13) 


dE _ 2E 
duj 3w ’ 

which is a differential equation for E(u}). The solution is 
E{uj) = (14) 


which is exactly the same asymptotic behavior as given 
by equation Eq. (11). In fact, the proportionality con¬ 
stant A can be found by approximating (r^) « Tq, where 
ro is the minimum of the effective potential. Then, 
A = is found, consistent with Eq. (11). 

Moreover, with the correct asymptote for w —> oo given 
by Eq. (9) and with a suitable two-parameter interpola¬ 
tion function we find a more complete differential equa¬ 
tion for general n and I, namely 


dE _2E n 
doj 3a; 3 


2ti ^ 2 


h 1 
”*"31 + auj ^ 




where a and /3 are fitting parameters. The frequency 
dependence of the last term is an empirical choice that 
interpolates between Taut’s solutions, for which (r^) can 
be calculated analytically. After fitting through the exact 
solutions, the above differential equation can be solved 
analytically to find the eigenenergy as function of fre¬ 
quency. Note that this provides us with a fundamen¬ 
tally different interpolation scheme than the one in [18]; 
whereas there the Hellman-Feynman theorem is used af¬ 
ter establishing E{lo) here it is used predictively. Another 
significant difference is that our method easily allows for 
interpolation of excited states (n 7 ^ 0 or Z 7 ^ 0 ). 


Interpolated eigenenergies 


-n = 0./ = 0 

- n = 0,/=1 

- n = 0,l = 2 

-^n = 0,l = 3 

-n = 0, / = 4 

- n=^,l = 0 

- n=^J=^ 

- n=1./ = 2 
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Figure 3: Eigenenergies for a number of states, with quantum 
numbers as indicated. The asymptotic behavior for a; —>■ 0 
has been subtracted for a clearer view on the results. 



The figure might indicate that for a; —0 the energy 
does not depend on the quantum number 1. With the 
full effective potential of Eq. (7), the harmonic approxi¬ 
mation Eq. (10), and expanding the resulting energy up 
to first order in /(I -h 1 ) we find 


nco nuj OLO 


(16) 


where is given by Eq. (11) and 


E^^\oj) 1 1 / 

hui 2 '^/^\U!o) 2) \UJo) 

(17) 

This result was found numerically in the context of har¬ 
monium [17], but here we derived it analytically, which 
allows for a scaling in terms of the Hartree frequency. 
Higher order corrections can also be found analytically. 
Now, since lini ~ *1’ energy does not depend 

any more on I in the limit of low frequencies, and is sim¬ 
ply given by Eq. (11). This was to be expected since 
we already had established that the importance of the 
Coulomb term in Eq. (7) is negligible, and therefore the 
centrifugal barrier term should be negligible as well. The 
limiting case in which the energy is independent of the 
quantum number /, indicates that states with different 
angular momenta are equally important. This is in agree¬ 
ment with the partial wave treatment of Coulomb scat¬ 
tering [ 21 ], which formally corresponds to a; = 0 . 


GENERAL STATES 


MANY-PARTICLE SYSTEMS 


Fitting Eq. (15) through the exact solutions found by 
Taut’s method we can determine the eigenenergies as 


A gas of degenerate electrons as mentioned in the in¬ 
troduction, is generally a multi-electron system. If we 
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assume only pairwise interactions between electrons, we 
can easily use the results from previous sections to find 
the ground state by populating the two-particle states 
given in Fig. 3. In the weakly-coupled regime this cor¬ 
responds to filling the eigenstates of a simple harmonic 
oscillator from below. In the strongly-coupled regime, for 
sufficiently small w, all pairs are in states with n = 0. In 
the limit of w —>■ 0, these states all become energetically 
degenerate which is of course a non-trivial situation for 
fermions. 

With the techniques shown in this paper, finding an¬ 
alytical expressions for the energies of eigenstates in a 
many-particle electron system should be possible. In the 
weakly-coupled regime the Coulomb interaction energy is 
much smaller than the harmonic oscillator energy separa¬ 
tion. A many-particle electron system behaves then as a 
weakly-interacting Fermi gas with the pairwise Coulomb 
interaction energy as a perturbation. Such an electron 
system can be quantum degenerate if its temperature is 
below the Fermi temperature of the harmonic oscillator. 
In the strongly-coupled regime where the Coulomb in¬ 
teraction energy dominates, perturbative methods seem 
at first sight not applicable. However, subtraction of the 
diverging energy scaling, Eq. (14), results in a har¬ 
monic oscillator problem with a factor y/S frequency and 
mixed angular momentum states, where the remaining 
Coulomb contribution is again small compared to the 
harmonic oscillator energy separation. This, combined 
with some of the results in the context of harmonium, 
suggests that also here a perturbative many-body ap¬ 
proach is possible. 

In fact, for N trapped electrons it can be shown [20] 
that the Hamiltonian is given by 


H = -- 


2 m« 


^ 1 


P=1 


N 

E 

p>q^l 


rrieUJ 

2N 


pq 


(18) 

N 

where R = ^pq = ^p-r^ and Vp = gf- with 

p=i ^ 

the position of the p-th electron. In the weakly-correlated 
regime, we neglect the Coulomb repulsion like we did for 
N = 2, so that the potential is purely harmonic, and 
therefore the relative energies are given by 


E 

huj 



(19) 


which indeed describes filling a simple harmonic oscilla¬ 
tor with electrons from below. An approximation for the 
energy in the strongly-correlated limit can easily be found 
if we assume that the coordinates Vpg in the Hamiltonian 
are independent. Then, making a harmonic approxima¬ 
tion around the minimum of each of the pair-potentials 
and by taking into account the multiplicity of the dou¬ 
ble sum, we can solve for this Hamiltonian in terms of 


^N{N —1) single-particle harmonic oscillator states with 
a total energy given by 




^ = ( — 
hio A \uJo 

^N{N- 

+ E 

i 

UJ \ 


-1/3 




2=1 

-1/3 


ni + 


+ B 


N 


^N(N-l) 

E 


( 20 ) 


( 21 ) 


The first term in Eq. (20) corresponds to a frozen-gas ap¬ 
proximation where the kinetic energy of all electrons is 
neglected and simply the minimum of each pair-potential 
is determined. The minimum of the potential results in 
an equilibrium position between the electrons, where the 
spatial positions form a crystal-like pattern. The cor¬ 
responding distance that minimizes the pair-potential is 
ro = {oj for any N. The electron gas in 

this regime is equivalent to a Wigner crystal [3, 19, 27] up 
to A^ = 4. The resulting coefficient = 3N'^^^{N—1)/A 
matches exactly the expressions found in Eq. (11) for 
N = 2 and in Refs. [19, 20] for A^ = 3 and A^ = 4. 
The second term in Eq. (20) corresponds to vibrational 
states of the crystal. Comparing the coefficient Hat = 
N{N — l)-\/3/4 to the results for A^ = 2 in Eq. (11) and 
N = 3,4 in Refs. [19, 20] we see that the relative differ¬ 
ence is 0%, 24% and 46% respectively. We interpret this 
difference as correlation energy due to the angular cor¬ 
relations between vibrational states of the crystal, which 
are correctly accounted for in Refs. [19, 20] but we ne¬ 
glected by assuming that all Ppq are independent. How¬ 
ever, far into the strongly correlated regime the result¬ 
ing energy corrections are relatively small compared to 
the total energy E. Thus we see that also for N elec¬ 
trons, harmonic approximations in both the weakly and 
strongly correlated regime indeed are suitable methods to 
find the eigenenergies and corresponding wave functions. 

The Hellman-Feynman interpolation method that we 
presented can also be generalized to many-particle elec¬ 
tron systems. However, in contrast to A^ = 2 where we 
used Taut’s exact solutions to determine the fitting pa¬ 
rameters a and /3 in Eq. (15), we would have to resort 
to numerical solutions as in Refs. [27, 28]. Eurthermore, 
to find the generalized form of Eq. (15), we can use the 
expressions Eq. (19) and Eq. (20), or if a higher accuracy 
is needed one could use the results from Refs. [19, 20]. 


CONCLUSION 

We described the problem of two electrons trapped in 
a harmonic potential for any value of the trapping fre¬ 
quency. Two different regimes can be identified: the 
strongly-coupled regime for small frequencies, and the 
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weakly-coupled regime for large frequencies. In both 
regimes we found a harmonic approximation of the effec¬ 
tive potential to correctly produce the eigenenergies of 
the electron pair for the ground state as well as any ex¬ 
cited state. The critical crossover trapping frequency cor¬ 
responds to the Hartree frequency. Using the Hellmann- 
Feynman theorem we were able to find semi-analytical 
solutions to the problem at any given frequency and for 
any given state. As shown, these concepts are readily 
generalized to many-particle systems. 

The trapping frequencies for current electron traps 
are several orders of magnitude below the crossover fre¬ 
quency. This implies that present day experiments are 
far into the strongly-coupled limit. Conversely, in the 
context of harmonium trapping frequencies are typically 
above the crossover frequency [16]. Interestingly, as we 
demonstrated in this letter, both these regimes can be 
described analytically, experimentally relevant to the ex¬ 
periments proposed on entangled electron pairs [6, 8]. 
Moreover, our methodology allows for a perturbative 
many-body calculation based on a paired-electron de¬ 
scription, as well as a feasible interpolation method. This 
allows us to explore the properties and phases of many- 
particle quantum-degenerate electron systems under har¬ 
monic confinement. 

This research was financially supported by the Foun¬ 
dation for Fundamental Research on Matter (FOM). 


* Electronic address: r.m.kroeze@student.tue.nl 

[1] APS, Accelerators and beams, tools of discovery and in¬ 
novation, Brochure (2010), URL http://www.aps.org/ 
units/dpb/news/edition4th.cfm. 

[2] O. Luiten, B. Claessens, S. Van der Geer, M. Reijnders, 
G. Taban, and E. Vredenbregt, International Journal of 
Modern Physics A 22, 3882 (2007). 

[3] E. Wigner, Physical Review 46, 1002 (1934). 

[4] W. Paul, Rev. Mod. Phys. 62, 531 (1990). 

[5] L. S. Brown and G. Gabrielse, Rev. Mod. Phys. 58, 233 
(1986). 


[6] L. Lamata, D. Porras, J. I. Girac, J. Goldman, and 
G. Gabrielse, Phys. Rev. A 81 , 022301 (2010). 

[7] J. L. Ghaloupka and D. D. Meyerhofer, Phys. Rev. Lett. 
83 , 4538 (1999). 

[8] P. Lougovski and H. Batelaan, Phys. Rev. A 84, 023417 

( 2011 ). 

[9] J. B. Balewski, A. T. Krupp, A. Gaj, D. Peter, H. P. 
Buchler, R. Low, S. Hofferberth, and T. Pfau, Nature 
502 , 664 (2013). 

[10] S. Schmid, A. Hrter, A. Frisch, S. Hoinka, and J. H. 
Denschlag, Review of Scientific Instruments 83 , 053108 
( 2012 ). 

[11] T. Busch, B.-G. Englert, K. Rzazewski, and M. Wilkens, 
Foundations of Physics 28 , 549 (1998). 

[12] J. Mentink and S. Kokkelmans, Physical Review A 79 , 
032709 (2009). 

[13] A. Samanta and S. K. Ghosh, Physical Review A 42, 
1178 (1990). 

[14] M. Taut, Physical Review A 48, 3561 (1993). 

[15] N. R. Kestner and O. Sinanolu, Physical Review 128 , 
2687 (1962). 

[16] J. Gioslowski and K. Pernal, The Journal of Chemical 
Physics 113 , 8434 (2000). 

[17] J. Karwowski and L. Cyrnek, Annalen der Physik 13 , 
181 (2004). 

[18] J. Gioslowski, The Journal of Chemical Physics 136, 
044109 (2012). 

[19] J. Gioslowski and K. Pernal, The Journal of Chemical 
Physics 125 , 064106 (2006). 

[20] J. Gioslowski and E. Grzebielucha, Phys. Rev. A 77 , 
032508 (2008). 

[21] J. J. Sakurai, Modem Quantum Mechanics (Addison- 
Wesley Publising Company, 1994). 

[22] B. H. Bransden and C. J. Joachain, Quantum mechanics 
(Pearson Education, 2000). 

[23] A. Turbiner, Physical Review A 50 , 5335 (1994). 

[24] A. Turbiner, Communications in Mathematical Physics 
118 , 467 (1988). 

[25] H. Hellmann, Leipzig: Franz Deuticke (1937). 

[26] R. P. Feynman, Physical Review 56, 340 (1939). 

[27] J. Gioslowski, K. Strasburger, and E. Matito, The Jour¬ 
nal of Chemical Physics 136 , 194112 (2012). 

[28] J. Gioslowski, K. Strasburger, and E. Matito, The Jour¬ 
nal of Chemical Physics 141 , 044128 (2014). 



